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Recent observation of the power spectrum of Cosmic Microwave Background (CMB) Radiation has exhibited that 
the flat cosmology is most likely. This suggests too large universal baryon-density parameter Q.\ > h 2 « 0.022 ~ 0.030 
to accept a theoretical prediction, Q,bh 2 < 0.017, in the homogeneous Big-Bang model for primordial nucleosynthesis. 
Theoretical upper limit arises from the sever constraints on the primordial 7 Li abundance. We propose two cosmo- 
logical models in order to resolve the descrepancy; lepton asymmetric Big-Bang nucleosynthesis model, and baryon 
inhomogeneous Big-Bang nucleosynthesis model. In these cosmological models the nuclear processes are similar to 
those of the r-process nucleosynthesis in gravitational collapse supernova explosions. 

Massive stars > IOMq culminate their evolution by supernova explosions which are presumed to be the most 
viable candidate site for the r-process nucleosynthesis. Even in the nucleosynthesis of heavy elements, initial entropy 
and density at the surface of proto-neutron stars are so high that nuclear statistical equilibrium favors production of 
abundant light nuclei. In such explosive circumstances many neutron-rich radioactive nuclei of light-to-intermediate 
mass as well as heavy mass nuclei play the significant roles. 
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I. Big-Bang Cosmology 

Recent progress in cosmological deep survey has clarified 
progressively the origin and distribution of matter and evolu- 
tion of Galaxies in the Universe. 

The origin of the light elements among them has been a 
topic of broad interest for its significance in constraining the 
dark matter component in the Universe and also in seeking for 
the cosmological model which best fits the recent data of cos- 
mic microwave background (CMB) fluctuations. This paper 
is concerned with neutrinos during Big-Bang nucleosynthesis 
(BBN). In particular, we consider new insights into the pos- 
sible role which degenerate neutrinos may have played in the 
early Universe 1 ' . 

There is no observational reason to insist that the universal 
lepton number is zero. It is possible, for example, for the in- 
dividual lepton numbers to be large compared to the baryon 
number of the Universe, while the net total lepton number is 
small L ~ B. It has been proposed recently 2 ' that models 
based upon the Affleck-Dine scenario of baryogenesis might 
generate naturally lepton number asymmetry which is seven to 
ten orders of magnitude larger than the baryon number asym- 
metry. Neutrinos with large lepton asymmetry and masses 
~ 0.07 eV might even explain the existence of cosmic rays 
with energies in excess of the Greisen-Zatsepin-Kuzmin cut- 
off. 3 ' It is, therefore, important for both particle physics and 
cosmology to carefully scrutinize the limits which cosmology 
places on the allowed range of both the lepton and baryon 
asymmetries. 

1. Lepton Asymmetric Big-Bang Model 

Although lepton asymmetric BBN has been studied in 
many papers 4 ' (and references therein), there are several dif- 
ferences in the present work: For one , we have included finite 



temperature corrections to the mass of the electron and pho- 
ton. 5 ' Another is that we have calculated the neutrino annihi- 
lation rate in the cosmic comoving frame, in which the M0ller 
velocity instead of the relative velocity is to be used for the in- 
tegration of the collision term in the Boltzmann equations. 6, 7) 

Neutrinos and anti-neutrinos drop out of thermal equilib- 
rium with the background thermal plasma when the weak re- 
action rate becomes slower than the universal expansion rate. 
If the neutrinos decouple early, they are not heated as the par- 
ticle degrees of freedom change. Hence, the ratio of the neu- 
trino to photon temperatures, T„/T 7 , is reduced. The biggest 
drop in temperature for all three neutrino flavors occurs for 
£ v ~ 10. This corresponds to a decoupling temperature above 
the cosmic QCD phase transition. £„ is the neutrino degener- 
acy parameter defined by = n v /T v , where \i v is the chem- 
ical potential and T„ is the neutrino temperature. Finite 
leads to a lepton asymmetric (L ^ 0) Universe. 

Non-zero lepton numbers affect nucleosynthesis in two 
ways. First, neutrino degeneracy increases the expansion rate. 
This increases the 4 He production. Secondly, the equilibrium 
n/p ratio is affected by the electron neutrino chemical poten- 
tial, n/p = exp{-(AM/T n ) - £„ e }, where AM is the 
neutron-proton mass difference and T„ Hp is the freeze-out 
temperature for the relevant weak reactions. This effect ei- 
ther increases or decreases 4 He production, depending upon 
the sign of £ Uc . 

A third effect emphasized in this paper is that T v jT 1 can 
be reduced if the neutrinos decouple early. This lower tem- 
perature reduces the energy density of neutrinos during BBN, 
and slows the expansion of the Universe. This decreases 4 He 
production. 

Figure 1 highlights the main result of this study, where we 
take = £„ T . For low Slb^o models, only the usual low 
values for £„ e and T are allowed. Between h^ Q w 0.188 
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Neutrino degeneracy can even allow baryonic densities up to 
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Fig. 1 Allowed values of £„ e and T for which the constraints 
from light element abundances are satisfied for values of 



b'«5() 



: 0.075, 0.1, 0.2, 0.3 and 1.0 as indicated. Note that 



^bfrJo = 4 x Q b h 2 , and h = hi 



2. Cosmic Microwave Background 

Several recent works 8-10 ' have shown that neutrino degen- 
eracy can dramatically alter the power spectrum of the CMB. 
However, only small degeneracy parameters with the standard 
relic neutrino temperatures have been utilized. Here, we have 
calculated the CMB power spectrum to investigate effects of 
a diminished relic neutrino temperature. 

The solid line on Fig. 2 shows a £7a = 0.4 model for 
which n = 0.78, where n is the power index of primordial 
fluctuations. This fit is marginally consistent with the data at 
a level of 5.2a. The dotted line shows the matter dominated 
Qjs. = best fit model with n — 0.83 which is consistent with 
the data at the level of 3<r. The main differences in the fits 
between the large degeneracy models and our adopted bench- 
mark model are that the first peak is shifted to slightly higher I 
value and the second peak is suppressed. One can clearly see 
that the suppression of the second acoustic peak is consistent 
with our derived neutrino-degenerate models. In particular, 
the MAXIMA- 1 results are in very good agreement with the 
predictions of our neutrino-degenerate cosmological models. 
It is clear that these new data sets substantially improve the 
goodness of fit for the neutrino-degenerate models. 9 ' More- 
over, both data sets seem to require an increase in the baryonic 
contribution to the closure density as allowed in our neutrino- 
degenerate models. 1 ' 13 ' 

3. Baryon Inhomogeneous Big-Bang Model 

The biggest advantage of the baryon inhomogeneous Big- 
Bang nucleosynthesis model 14-16 ' is to allow larger Q\ } h 2 < 
0.05, which well covers the constraint from recent CMB data 
n b h 2 w 0.022 ~ 0.030, still satisfying the light element 
abundance constraints. 
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Fig. 2 CMB power spectrum from BOOMERANG 10 (squares) and 
MAXIMA- 1 12 ' (circles) binned data compared with calcu- 
lated Q — 1 models. 



Let us consider what kind of observational signature is ex- 
pected in this model. Nuclear reaction flow stops at the A 
= 7 nuclear systems in the homogeneous Big-Bang nucle- 
osynthesis model because of the instability of 8 Be. In the 
baryon inhomogeneous model, however, the nucleosynthe- 
sis occurs in an environment of proton-neutron segregated 
inhomogeneous distribution. Therefore, the radioactive nu- 
clear reactions play the significant roles in the production 
of intermediate-to-heavy mass elements via unstable nuclei 
8 Li(838 ms), 9 Li(178.3 ms), 7 Be(53.29 d), 10 Be(1.51 x 10 6 y), 
8 B(770 ms), 12 B(20.20 ms), n C(20.385 m), 14 C(5730 y), 
15 C(2.449 s), 13 N(9.965 m), 16 N(7.13 s), 14 O(70.606 s), etc. 



4 He( 3 H, 7) 7 Li(n, 7 ) 8 Li(a, n) n B(n, 7 ) 12 B(/3 v) 
12 C( n , 7 ) 13 C(n, 7 ) 14 C..., 



(1) 



7 Li(n, 7) 8 Li(n, 7) 9 Li(/5 ^) 9 Be(n, 7 ) 10 Be(n, 7) (2) 
n Bc(/3 v) n B..., 



7 Li( 3 H,n) 9 Be( 3 H,n) n B. 



(3) 



The two reaction chains (1) and (2) play the key roles in the 
production of heavy neutron-rich isotopes 14 ' 16 ' in the neutron- 
rich zones. Since 7 Li is the heaviest element to be created 



in the homogeneous Big-Bang nucleosynthesis model, all nu- 
clear reactions for the production of heavier elements have 
ever been ignored in the previous calculations. We however 
found that the reaction chain (3) is extremely important for 
the production of 9 Be in the baryon inhomogeneous Big-Bang 
nucleosynthesis models. 17 ~ 19) With these reactions being in- 
cluded in the network of the baryon inhomogeneous Big-Bang 
models, the 9 Be abundance increases by three orders of mag- 
nitude to N(Be)/N(H) « 1CT 14 which approaches the current 
observational level. 

II. Supernova Explosion 

Stars with various masses provide a variety of production 
sites for intermediate-to-heavy mass elements. Very massive 
stars > 10M© culminate their evolution by supernova (SN) 
explosions which are presumed to be most viable candidate 
for the still unknown astrophysical site of r-process nucle- 
osynthesis. We discuss in this section the neutrino-driven 
winds from Type II SN explosion of very massive stars. Al- 
though there is still a room for the prompt explosion 20 ' to ac- 
count for one part of the r-process nucleosynthesis, we con- 
centrate on the gravitaional core-collpase Type II SNe here. 

Even in the nucleosynthesis of heavy elements, initial en- 
tropy and density at the surface of proto-neutron stars are so 
high that nuclear statistical equilibrium (NSE) favors produc- 
tion of abundant light nuclei. In such explosive circumstances 
of so called hot-bubble scenario, not only heavy neutron rich 
nuclei but light unstable nuclei play a significant role. 

The study of the origin of r-process elements is also critical 
in cosmology. It is a potentially serious problem that the cos- 
mic age of the expanding Universe derived from cosmological 
parameters may be shorter than the age of the oldest globular 
clusters. Since both age estimates are subject to the uncer- 
tain cosmological distance scale, an independent method has 
long been needed. Thorium, which is a typical r-process ele- 
ment and has half-life of 14 Gyr, has recently been detected 
along with other elements in very metal-deficient stars. If we 
model the r-process nucleosynthesis in these first-generation 
stars, thorium can be used as a cosmochronometer completely 
independent of the uncertain cosmological distance scale. 

1. Neutrino-Driven Winds in Type-II Supernovae 

Recent measurements using high-dispersion spectrographs 
with large Telescopes or the Hubble Space Telescope have 
made it possible to detect minute amounts of heavy elements 
in faint metal-deficient ([Fe/H] < -2) stars. 21 ' The discovery 
of r-process elements in these stars has shown that the relative 
abundance pattern for the mass region 120 < A is surprisingly 
similar to the solar system r-process abundance independent 
of the metallicity of the star. Here metallicity is defined by 
[Fe/H] = log[N(Fe)/N(H)] - log[N(Fe)/N(H)] s . It obeys the 
approximate relation t/10 10 yr ~ 10^ Fe/>H ^. The observed sim- 
ilarity strongly suggests that the r-process occurs in a single 
environment which is independent of progenitor metallicity. 
Massive stars with 10M Q < M have a short life ~ 10 7 yr 
and eventually end up as violent supernova explosions, eject- 
ing material into the intersteller medium early on quickly from 



the history of the Galaxy. However, the iron shell in SNe is 
excluded from being the r-process site because of the observed 
metallicity independence. 

Hot neutron stars just born in the gravitational core col- 
lapse SNell release most of their energy as neutrinos during 
the Kelvin-Helmholtz cooling phase. An intense flux of neu- 
trinos heat the material near the neutron star surface and drive 
matter outflow (neutrino-driven winds). The entropy in these 
winds is so high that the NSE favors a plasma which consists 
of mainly free nucleons and alpha particles rather than com- 
posite nuclei like iron. The equilibrium lepton fraction, Y e , is 
determined by a delicate balance between v e +n — > p+e~ and 
v e +P — > n + e + , which overcomes the difference of chemical 
potential between n and p, to reach Y e ~ 0.45. R-process nu- 
cleosynthesis occurs because there are plenty of free neutrons 
at high temperature. This is possible only if seed elements 
are produced in the correct neutron-to-seed ratio before and 
during the r-process. 

Although Woosley et al. 22) demonstrated a profound pos- 
sibility that the r-process could occur in these winds, several 
difficulties were subsequently identified. First, independent 
non relativistic numerical supernova models 23 ' have difficulty 
producing the required entropy in the bubble S/k ~ 400. Rel- 
ativistic effects may not be enough to increase the entropy 
dramatically. 24 ~ 26) Second, even should the entropy be high 
enough, the effects of neutrino absorption v e +n — > p+e~ and 
v e + A(Z, N) -> A(Z+l,N-l) + e~ may decrease the neu- 
tron fraction during the nucleosynthesis process. As a result, 
a deficiency of free neutrons could prohibit the r-process. 27) 

In order to resolve these difficulties, we have studied 26 ' 28) 
neutrino-driven winds in a Schwarzschild geometry under the 
reasonable assumption of spherical steady-state flow. The pa- 
rameters in the wind models are the mass of neutron star, M, 
and the neutrino luminosity, L v . The entropy per baryon, 
S/k, in the asymptotic regime and the expansion dynamic time 
scale, Tdyn, which is defined as the duration time of the a- 
process when the temprature drops from T « 0.5 MeV to 0.5/e 
MeV, are calculated from the solution of hydrodynamic equa- 
tions. Then, we carried out r-process nucleosynthesis calcu- 
lations in our wind model. We found 26 ' that the general rela- 
tivistic effects make Td yn much shorter, although the entropy 
increases by about 40 % from the Newtonian value of S/k ~ 
90. By simulating many supernova explosions, we have found 
some interesting conditions which lead to successful r-process 
nucleosynthesis, as to be discussed in the following sections. 

2. R-process Nucleosynthesis 

Previous r-process calculations 22 ' 30) had complexity that 
the seed abundance distribution was first calculated by using 
smaller network for light-to-intermediate mass elements, and 
then the result was connected further to another r-process net- 
work in a different set of the computing run. For this reasaon 
it was less transparent to interpret the whole nucleosynthesis 
process. This inconvenience happened because it was numer- 
ically too heavy to run both a-process and r-process in a sin- 
gle network code for too huge number of reaction couplings 
among ~ 3000 isotopes. Our nucleosynthesis calculation 26 ' 28 ' 
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Fig. 3 R-process abundance 26 ' (solid line) as a function of atomic 
mass number A compared with the solar system r-process 
abundance (filled circles) from Kappeler, Beer, & Wisshak. 29 ' 
The neutrino-driven wind model used is for L v — 10 52 ergs/s 
and M — 2M Q . The solar system r-process abundance is 
shown in arbitrary unit. 
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Fig. 4 The same as those in Figure 3, but for the neutrino-driven 
wind model of L v — 5 x 10 5 ergs/s. Solid line respre- 
sents the result by using the Woosley & Hoffman rate 31 ' of 
the 4 He(an, 7) 9 Be reaction, and long-dashed line for the rate 
multiplied by factor 2, as suggested by the recent experiment 
of Utsunomiya et al. 32) 



is completely free from this complexity because we exploited 
fully implicit single network code which is applied to a se- 
quence of the whole processes of NSE - a-process - r-process. 



Let us remind the readers that there were at least three dif- 
ficulties in the previous theoretical studies of the r-process. 
The first difficulty among them is that an ideal, high entropy 
in the bubble S/k ~ 400 22 ' is hard to be achieved in the other 
simulations. 23-26 ' 

The key to resolve this difficulty is found with the short 
dynamic time scale Td yn <~ 10 ms in our models 26 ' 28) of 
the neutrino-driven winds. As the initial nuclear composi- 
tion of the relativistic plasma consists of neutrons and pro- 
tons, the a-burning begins when the plasma temperature cools 
below T ~ 0.5 MeV. The 4 He(aa,7) 12 C reaction is too 
slow at this temperature, and alternative nuclear reaction path 
4 He(cm, 7) 9 Be(a, n) 12 C triggers explosive a-burning to pro- 
duce seed elements with A ~ 100. 3I) Therefore, the time scale 
for nuclear reactions is regulated by the 4 He(an, 7) 9 Be. It is 
given by tat = {plY^Y n \{aan -^ 9 Be)) . If the neutrino- 
driven winds fulfill the condition Td yn < tjv, then fewer 
seed nuclei are produced during the a-process with plenty 
of free neutrons left over when the r-process begins at T ~ 
0.2 MeV. The high neutron-to-seed ratio, n/s ~ 100, leads 
to appreciable production of r-process elements, even for low 
entropy S/k ~ 130, producing both the 2nd (A ~ 130) and 
3rd (A ~ 195) abundance peaks and the hill of rare-earth ele- 
ments (A ~ 165) (Fig. 3). 

The three body nuclear reaction cross section for 
4 He(cm, 7) 9 Be is one of the poorly determined nuclear data 
which may alter the r-process nucleosynthesis yields. The 
inverse process has recently been studied experimentally by 
Utsunomiya et al., 32) and photodisintegration cross section of 



9 Be has been measured with better precision than those of the 
previous experiments. Applying the principle of the detailed 
balance to this process, one can estimate the cross section for 
4 He(cm, 7) 9 Be. They found that the thermonuclear reaction 
rate is almost twice as big as that of Woosley and Hoffman 31 ' 
but in resonable agreement with recent compilation of An- 
gulo et al. 33 ' However, there still remain several questions on 
the consistency between their result and electron-scattering 
experiments, on the contribution from the narrow resonance 
,r = 5/2" (2.429 MeV), etc. It is also a theoretical chal- 
lenge to understand the reaction mechanism and the resonance 
structure because two different channels, 8 Be + n and 5 He + 
a, contribute to this process. 

Therefore, we show two calculated results in Fig. 4: The 
solid line displays the result obtained by using the Woosley 
and Hoffman cross section, 31 ' assuming a 8 Be + n structure 
for 9 Be. We also calculated the r-process by multiplying this 
cross section by factor of 2 (long-dashed line). This makes a 
drastic change in the r-process yields in the 3rd (A ~ 195) 
abundance peak. More theoretical and experimental studies 
of the 4 He(cm, 7) 9 Be reaction are highly desired. 

3. Neutrino-Nucleus Interactions 

Neutrino interactions with nucleons and nuclei take the key 
to resolve the second difficulty which was pointed out in sect. 
1. The difficulty is that the effects of neutrino absorptions 
v e + n -^p + e- and v e + A(Z, N) -> A(Z + 1, N- 1) + e~ 
during the a-process may induce the deficiency of free neu- 
trons and break down the r-process conditions. 27 ' These 
two types of neutrino interactions control most sensitively 
the electron fraction and the neutron fraction, as well, in a 
neutron-rich environment. In order to resolve this difficulty, 
we have updated the electron-type neutrino capture rates for 
all nuclei and electron-type anti-neutrino capture rate for free 
proton. 34 ' 35 ' 



The new r-process calculation proves to be almost invari- 
ant. One can understand this robustness of the succesful r- 
process in the following way: The specific collision time for 
neutrino-nucleus interactions is given by 

(4) 

where Li. 51 is the individual neutrino or antineutrino lumi- 
nosity in units of 10 51 ergs/s, =< Ef > / < E t > in 
MeV (i = u e , D e , etc.), and (<j v ) is the averaged cross 
section over neutrino energy spectrum. At the a-burning 
site of r « 100 km for L„ >51 w 10, e Ve = 12 MeV, and 
(o-„) w 10~ 41 cm 2 , T„ e (r=100 km) turns out to be w 240 ms. 
This collision time is larger than the expansion dynamic time 
scale; Td yn ~ 10 ms <C 7v e (r=100 km) w 240 ms. Because 
there is not enough time for v e 's to interact with n's in such 
rapidly expanding neutrino-driven wind, the neutron fraction 
is insensitive to the neutrino absorptions. 

One might wonder if our dynamic time scale ~ 10 ms is too 
short for the wind to be heated by neutrinos. Careful compari- 
son between proper expansion time and specific collision time 
for the neutrino heating is needed in order to answer this ques- 
tion. Otsuki et al. 26) have found that the supernova neutrinos 
transfer their kinetic energy to the wind most effectively just 
above the neutron star surface at 10km < r < 20km. There- 
fore, one should refer the duration time for the wind to reach 
the a-burning site, Th oa t, rather than Td yn - One can estimate 
this heating time scale 

r r f d r 

Thcat = / , (5) 

where u is the fluid velocity of the wind. By setting the radius 
of neutron star surface r, = 10 km and r/ = 100 km, we get 
Thcat ~ 30 ms. The collision time r„ is given by Eq. (4) by 
setting L„ >B1 » 10, e„ = (e„. + e p J/2 = (12 + 22)/2 = 17 
MeV, r «1 km, and (er„) « 10~ 41 cto 2 . Let us compare Th ea t 
and t v to one another: t v w 3.4ms <C Thoat ~ 30ms. We 
can thus conclude that there is enough time for the expanding 
wind to be heated by neutrinos even with short dynamic time 
scale for the a-process, Td yn ~ 10 ms. 

4. Roles of Light Neutron-Rich Nuclei 

The r-process is thought to proceed after the pile up of seed 
nuclei produced in the a-process at higher temperatures Tg w 
5 ~ 2.5. Since charged-particle reactions, which reassem- 
ble nucleons into a-particles and a-particles into heavier nu- 
clei (i.e. a-process), are faster than the neutron-capture flow 
which is regulated by beta-decays, the light-mass neutron-rich 
nuclei were presumed to be unimportant. 

However, Terasawa et al. 36) have recently found that even 
light neutron-rich nuclei progressively play the significant 
roles in the production of seed nuclei. Nuclear reaction net- 
work used in the previous studies 22 ' 30) includes only limited 
number of light unstable nuclei, 3 H, 7 Be, 8 < 9 B, n < 14 C, 13 N, 
15q i8,20p^ 23,24 Ne ^ and SQ on We therefore need to extend 

the network code so that it covers all radioactive nuclei to the 



neutron-drip line. We take the rates of charged particle reac- 
tions from those used in the Big-Bang nucleosynthesis calcu- 
lations 16 - 18 - 19) and the NACRE compilation. 33 ' 

Let us briefly discuss preliminary result of the r-process 
calculation, using the extended reaction network. 36) At early 
epoch of the wind expansion, t < a few dozens ms, both 
temperature and density are so high that the charged par- 
ticles interact with one aother to proceed nucleosynthesis 
around the /3-stability line in the light-mass region A < 20. 
There are plenty of protons and a-particles as well as neu- 
trons at this epoch, and the main reaction flow is triggered by 
4 He(<m, 7 ) 9 Be: 22 - 31 > 

4 He(cm, 7) 9 Be(a, n) 12 C(n, 7) 13 C(a, n) (6) 
16 O(n, 7 ) 17 O(a,n) 20 Ne or 16 0(a, 7 ) 20 Ne... 

However, at relatively later epoch even after the a-rich freeze 
out, a new reaction path 36) 

3 H(a, 7 ) 7 Li(n, 7) 8 Li(a, n) n B(n, 7 ) 12 B(n, 7) (7) 
13 B(n, 7 ) 14 B(n, 7) 15 B( C -^) 15 C... 

also takes some appreciable flux of baryon number to con- 
tinuously supply the seed nuclei. The classical r-process like 
flow, (n,7) followed by beta decay, has already started from 
light nuclei. This is a very different result from the previous 
picture that the r-process starts from only intermediate-mass 
seed nuclei A w 100. 

Since we do not have much information of (2n, 7) reac- 
tions, we did not include 6 He, 8 He, n Li, 14 Be, 17 > 19 Be, 22 C, 
etc. The yields of even the most neutron-rich isotopes were 
found to be abundant in this calculation, 36 ' and we plan to 
study the possible role of the (2n,7) reactions. There are sev- 
eral branching points between (n,7) and (a,n) reactions. They 
are at 18 C, 24 0, 36 Mg, etc. Exerimetal studies to measure 
these reaction cross sections are highly desirable. 

III. Quest for Nuclear Data and Astrophysics Data 

Our result in Fig. 3 reproduces fairely well the observed 
second abundance peak (A w 130), the hill of rare-earth el- 
ements (A w 165), and the third peak (A w 195). However, 
there are several defects, too. The first defect is a shift of 
the third abundance peak around Aw 195 by a couple of 
mass units. This is a common feature found in the previous r- 
process calculations, 22 ' 26 ' 30) too. These elements are the beta- 
decay products of extremely neutron-rich unstable nuclei on 
the neutron magic N = 126. Peak position depends on the 
timing of freezeout of the r-process. Therefore, a particular 
combination of environmental evolution of neutron-number 
density, N n , and temperature, Tg, as well as the expansion 
dynamic time scale, Td yn , might match the freeze out so that 
it results in the right position of abundance peak. 36) 

The second defect is the deficiency of abundance right 
above or below the peak elements, i.e. at A w 90, 120, 150, 
190, and 210. These deficiencies seem related to yet unseen 
effects of deformation or strucure change of unstable nuclei 
surrounding the neutron magic numbers N = 50, 82, and 126. 



Further extensive theoretical studies and observational chal- 
lenge to determine the masses, lives, and beta Q-values of 
these nuclei are highly desired. 

The third defect is the underproduction of actinoid ele- 
ments, Th-U-Pu (A = 230 ~ 240), by more than one order of 
magnitude. The observed high abundance level of these nuclei 
might suggest an existence of a new magic number around N 
= 150 - 160: Xenon 129 Xe 75 and platinum 195 Ptn 7 are the 
typical r-process elements on the second and third abundance 
peaks, which are decay products from extremely neutron-rich 
unstable nuclei with neutron magic numbers N = 82 and 126, 
respectively. From these observations we estimate that the 
waiting point nucleus is located by shifting AN w 7 or 9 units 
from the peak element. Applying the same shift AN w 7 or 
9 to 232 Thi42 and 238 Ui46, we could assume a new magic 
number around N = 150 ~ 160 which may lead to the fourth 
abundance peak at A = 230 ~ 240. Actually, in the very light 
nuclear systems, a new magic number N = 16 was found 37 ' 
in careful experimental studies of the neutron separation ener- 
gies and interaction cross sections of extremely neutron-rich 
nuclei. Since these possibilities were not taken into account in 
the present and previous calculations, the deficiency of acti- 
noids might be improved by modernizing nuclear mass for- 
mula including such effects. Another possibility is to make 
actinoid elements in neutron star mergers or the mergers of 
the neutron star and black hole binaries which have extremely 
small lepton fraction, Y e < 0.2. 38) However, these processes 
do not virtually produce any intermediate-mass nuclei includ- 
ing iron, which contradicts with the fact that the observed iron 
abundance is proportional to the r-process and actinoid ele- 
ments over the entire history of Galactic evolution ~ 10 10 yr. 

Since we discuss only material ejected from the proto- 
neutron star behind the shock, it does not make any serious 
problem to see the underproduction in mass region A < 90. 
Most of these intermediat-mass nuclei are ejected from the 
exploded outer shells in supernovae. 
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